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ABSTRACT 

I present the results of a search for intergalactic hydrogen clouds in voids. Clouds are detected by 
their HI Lya absorption lines in the HST spectra of low-redshift AGN. The parameter with which the 
environments of clouds are characterized is the tidal field, for this places a lower limit on the cloud mass- 
density which is dynamically stable against disruption. The summing of the tidal fields along these lines 
of sight is managed by the use of galaxy redshift catalogs. The analytical methodology employed here has 
been designed to detect gas clouds whose expansion following reionization is restrained by dark matter 
perturbations. The end-products of the analysis of data are the cloud equivalent width distribution 
functions (EWDF) of catalogs formed by sorting clouds according to various tidal field upper, or lower 
limits. 

Cloud EWDFs are steep in voids (d\og{dN / 'dz) / 'rflog(W) =S 1.5 ± 0.2), but flatter in high tidal 

field zones (S ~ —0.5 ± 0.1). Most probable cloud Doppler parameters are ~ 30 kms -1 in voids and 
~ 60 kms" 1 in proximity to galaxies. In voids, the cumulative line density at low EW (W ~ 15 mA) is 
essentially equal to that of the mean EWDF, ~ 500 per unit redshift. The void filling factor is found to 
be 0.87 < f v < 0.94. The void EWDF is remarkably uniform over this volume with a possible tendency 
for more massive clouds to be in void centers. The size and nature of the void cloud population suggested 
by this study is completely unanticipated by the results of published 3-D simulations, which predict that 
most clouds are in filamentary structures around galaxy concentrations, and that very few observable 
absorbers would lie in voids. Strategies for modeling this population are briefly discussed. 

Subject headings: intergalactic medium — quasars:absorption lines - dark matter - galaxies :halos 



1. INTRODUCTION 

Suppose we divide the universe into two complementary 
parts; one galaxy-dominated, and the other what might be 
called space-dominated, or void. While we know the pri- 
mary constituents of galaxy dominated space (hereafter, 
CDS), little is known of voids. What lies in voids? The 
presence of Lya absorbers in low-redshift spectra raises 
the possibility that some may be in voids. What is the 
nature of Lya absorbing clouds, and how does their na- 
ture change as one looks from void to GDS? What might 
be the long-term physical interaction between voids and 
the GDS? What remnants of primordial perturbations re- 
main in voids to this day? What is the indicative total 
mass density in voids? How does this compare to predic- 
tions based on numerical simulations? These are questions 



whichl motivated the present investigation Tot ns take a 

look at the clouds of the local universe; let us see what 
they are doing, with particular attention to those which 
are isolated. 

1.1. Clouds Observed 

To begin the study of the effect of environment on 
the nature of intergalactic hydrogen clouds one needs the 
equivalent width (EW) of their absorption at the Lya res- 
onance (A Q =1215.67 A), and their line profiles, charac- 
terized by the Doppler parameter b (kms -1 ), needed to 
calculate the column density of HI from the equivalent 
width. In addition, we shall require a parameter that will 
discriminate between void, and GDS clouds (to be dis- 
cussed in §2). 

The cloud population is characteristically measured by 
its line density dAf /dz, where Af is the number of clouds 



with column density greater than some Nhi- The line 
density may vary with redshift. This evolution has been 
expressed as a power of (1 + z); dAf /dz — (dAf/dz) (l + zp 
( Sargent et al. 1980 , hereafter SYBT), where 7 is some- 
times known as the evolution parameter. Before the 
Hubble Space Telescope (HST) made observation of low- 
redshift Lya clouds possible, our ideas about the local pop- 
ulation of absorbers were based on extrapolations from 
the observations of highly redshifted absorption spectra. 
These studies showed a strong evolution of the line den- 
sity dAf /dz with timeQ. 

A series of papers from the HST Quasar Absorption 
Line Key Project provided information on clouds at red- 
shifts z < 1.5 using the Faint Object Spectrograph (FOS) 
jBahcall ct al. 199l| , |1992| , |1993| , |1996j ; |Wcymann ct al 
1998| ). These papers provided palpable proof that the evo- 



lution rate of Lya clouds at low redshift is substan tially 
lower than at higher redshift. Wcymann et al. (1998) show 
that for clouds with no associated metal-line absorbers, the 
data is consistent with 7 = at z < 1. 5. For HI clouds 
with associated metal lines, 7 is in the range of 1.24 to 
1.55, depending on sample selection criteria. These ab- 
sorbers have been interpreted as arising from galactic halos 
(SYBT). 

Weymann et al. (1998 1 ) demonstrated that at low red- 
shifts 7 declines with EW (see their Fig. 7) indicating 
that low-EW clouds are positively evolving (i.e., becom- 
ing more numerous) relative to larger absorbers. Hl-only 
clouds with large EW (^ 400 mA) were found to have evo- 

1 Lu et al. (1991 ); Kim et al (1 qqTjl 
±0.71 respectively, though ticchtold (] 



1 ±0.28. 



show 7 ~ 2.75±0.29 and 2.78 
tiechtold (1994J) found a slower evolution 
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lution parameters 7 ~ 0.6, while the smaller clouds with 
W = 350 to 100 mA have 7 ~ 0.5 to —0.5, respectively. 
The evolution parameter of very low EW clouds at low red- 
shift is still unknown, but t he large redshift frequency of 



low EW absorbers found in Pcnton, Stocke, & Shull 2000 



(hereafter PSSI) suggests that the correlation between 7 
and EW may continue as the EW approaches 10 mA. 

The HST has brought the study of Lya absorbers down 
to a redshift at which cloud environments might b e stud- 
ied. An ea rly study of the sightline to 3C 273 (Morris 
et al. 1993 ) pointed out the tendency of clouds to avoid 
clusters of galaxies, and commented on the existence of 
one larg e but very isolated c l oud ^. Subsequently, other 
groups (Lanzetta et al. 1995 , Tripp, Lu, fc Savage 1998 



Chen, Lanzetta, Webb, fc Barcons 1998 ) investigated the 
proximity of clouds to galaxies, finding an anticorrelation 
between cloud equivalent width and the projected radius 
to the nearest galax y. The same relationship was found to 
hold out to z fa 0.8 ( |Chen et al., 1998| , |200l| ). [Tripp et al " 



(1998(J found that the anticorrelation may extend to pro 
jected radii of ~ 2 Mpc from the nearest galaxy, suggesting 
that galaxies may affect cloud distributions at very large 
galactocentric distances. Tripp et al. (1998 ) claimed that 
clouds beyond this radius had a distribution of equivalent 
widths similar to those clustered around galaxies, but that 
"void" clouds appear to be generally smaller. While galac- 
tic impact parameters within ~ 600 kpc were found to 
reliably produce absorption systems, outside this radius, 
increasingly frequent 'misses' were recorded, implying a 
declining covering factor with increasing radius at a given 
sensitivity limit. 

The view that these clou ds are physically associated 
with galaxies is challenged by Dave et al. 1999 . They claim 
that the absorption systems are associated with uniformly 
dense filamentary structures seen in 3-D hydrodynamical 
simulations, rather than self-gravitating clouds physically 
associated with galaxies and groups of galaxies. Whether 
associated with galaxies or filaments, there does appear to 
be a correlation showing the effects of proximity to galax- 
ies. The large physical exte nt of the radial trend in EWs 
noted by Tripp et al. (1998 ), and the impression of an in- 
creasing cloud number density toward galaxies, suggests 
that the surrounding void edges are the source of clouds 
which are falling into the gravitational potential wells of 
galaxies and groups of galaxies. A very similar distribution 
of clouds is detected out to redshifts z ~ 0.89 (Chen et al. 



1998 |, |200l| ), suggesting that such infall is a long-term fea- 
ture of the universe. Other lines of evidence favoring a 
substantial accretion rate of centrally condensed clouds 
are to be found in the roughly st eady, or increasing star 



formation rate in the Galaxy disk ( Maciel fc Rocha-Pinto, 
19980 the large dispersion of metallicities of disk stars as 



a fun ction of the epoch of their formation (Edvardsson 
et al.J 1993; Carraro et al., 1998), photometric asymme- 



tries and ot her signs of accretion of low-mctallicity gas in 
field spirals (Zaritsky, 1995; Zaritsky fc Rix, 1997; Hayncs 



200l|] the predominance of negative barocentric veloci 
ties for high velocity clouds (HVC) i n the Local group 
(Blitz et al., 1999; Wakker et al., 199E), and the existence 



2 The cloud in question is the highest wavelength cloud (point with 
error-bar) in Fig. 6a, in the sightline to 3c273, with an equivalent 
width of 297 mA 



of a "compact" subclass of infalling (~ 100 barocentric 
kms' 1 ) HVCs with cold , compact cores and ionized halos 
flBraun fc Burton, 2000[). 



1.2. Hydro Simulations 

The current view of the nature of the structure and con- 
tent of the universe at low redshift is strongly affected 
by results of numerical cosmological simulations. Vari- 
ous 3-dimensional simulations have been used to model 
the higher redshift Lya cloud p o pulation {e.g., pTernquist 
et al. 1996|; iMficket et al., 1996] ; pave et al., 199^ ; |Ccn g 



low redshift ([Riediger, Pctitjean, fc Muckct, 1998 j_ The- 



Ostrikcr. 19991 ) , and some of them have been extended to 



uns, Leonard, fc Efstathiou, 1998| ; pave et al., 1999 ) 
These simulations suggest that the clouds with columns 
log Nhi 14 are produced in the filaments of roughly ho- 
mogeneous gas, some of which may still be in partia l par- 
ticipation with t he Hubble flow (|Ricdiger et al., 199*8 ; The 
uns et al., 19981 Pave et al., 1999 ) , whereas the higher col- 



umn density clouds arc found in close proximity to galaxy 
halos. 

Very few of the simulations concern themselves with 
voids per se, so some effort mu st be taken to extract pre- 
dictions. iRiediger et al. (1998) found it helpful to sepa- 



rate clouds into two groups - shocked P s and unshocked 
P u . Shocked clouds result from convergent gas flows. It 
is thought that the unshocked clouds lie outside of the fil- 
amentary structures, and are expanding with the Hubble 
flow so that their column densities would decline steeply 
with declining redshift (Nhi oc (1 + z) 5 ) at higher red- 
shifts where J(z) is fairly constant (Riediger et al 



1998 



Dave et al., 1999] ) , so that this c omponent produces neg- 



Cen fc Ostriker (1999[ ) 
baryons are in the fil- 



ligible absorption for z < 3.0 
argue that from 1/2 to 2/3 of al 
amentary structures constituting a warm to hot medium 
(10 5 K to 10 7 K) with a volume filling factor ~ 0.1. T his 
is the population referred to by |Riediger et al. (1998 ) as 
"shocked" , and by Dave et al. (199~9~ ) as producing columns 
log Nhi <, 14 ( cm ) at low redshift. Void clouds are not 
expected to produce any significant Lya absorp tion (Dave, 
priv. comm . 2001). However, a simulation by Cen fc Os] 



trikcr (1999) finds that about 26% of baryons are located 
in a "warm" medium (T < 10 5 K) with a filling factor of 
0.9 - this would be the "unshocked" population P u whose 
contribution to the line density faded rapidly at redshifts 
less than 3. Thus, these simulations suggest that the un- 
shocked region is what we would call "void" , and has a 
current density of pjj w 0.26 of the average, and that 
filaments, with filling factor 0.1, have a density ps J> 7 
times the average, so that the void matter density would 
be less than 4% of that in the filaments. The issue of the 



void m atter density is addressed directly by Einasto et al i 
|(1994| ) through simulations of CDM cosmogony, showing 
that matter flows out of low density environments and into 
high density areas, resulting in an es timated 15% of matter 
being located in present-day voids. El-Ad & Piran (1997) 
estimate a void matter density of 10% of the mean based 
on observation of the galaxy content of nearby voids. 

Given the above, it seems fair to conclude that the pre- 
dictions of simulations are for isolated (i.e., void) clouds 
to be evolving rapidly, and for the line density in voids, 
at EWs detectable with current technonogy, to be quite 
low. The overwhelming majority of absorption systems 
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are expected to be produced by gas associated with the 
filamentary structures having warm to hot temperatur es 
and a filament filling factor // ^0.1 ( Dave et al., 1999 ). 

1.3. Structure Formation 

The conclusions from more theoretically-based stud- 
ies of structure formation are somewhat different. In 



semi-analyt i c simulations of very high rcdshifts ( Aba 



et al. 1998 Bromm, Coppi, fc Larson 1999 ; Nakamura 



& Umemura 2001), molecular hydrogen, formed shortly 
after recombination, provides cooling for low-mass pertur- 
batio ns which results in the formation of Population III 



stars |at z ~ 20 in halos with halo velocities as low at 
~ 1.5 kms -1 suggesting that even very small halos with 
masses in the range 10 5 to 10 9 M Q , would have formed 
stars. However, this field currently suffers from a "missing 
halos" problem; there should be many more dwarf galax- 
ies within the Galaxy's halo under the favored cold dark 
matter cosmogony (CDM) than are in fact observed. Es- 
timates of the magnitude of the discrepancy between pre- 
dicted and observed low -mass galaxies rang e from a few 
( Kauffmann ct al., 1993 ) to a factor of ~ 50 (Klypin et al 



1999| Moore et al. 1999a). It appears that the discrepancy 



manifests itself increasingly below halo velocities of ~ 40 
kms . The deficit in numbers between predicted low- 
mass mini-galaxies, and that suggested by observations of 
Local Group galaxies implies that the overwhelming ma- 
jority of small halos (v c <, 25 km s -1 ), for whatever reason, 
did not in fact generally form stars. Therefore, these halos 
could very well be present in the universe today as clouds. 
But while as mini-galaxies they could survive the tidal 
fields around large galaxies, as clouds the ram pressure 
of dissipated gas around galaxies and groups of galaxies 
would cause the stripping of halos of their baryons, ren- 
dering them invisible. Therefore, the cloud density may be 
less than the numbers of mini-galaxies predicted to have 
survived to the present with identical halo velocities due 
to their disruption in dense environments. Their remnants 
may remain as isolated (Lya clouds) , or they may be cur- 
rently falling into galaxy halos an d ide ntified as HVCs as 
suggested by Klypin ct al. (1999[) and |Grebcl. Braun, & 
Burton 20001- Yet it is claimed flKlypin ct al., 1999|) that 



their numbers are not inconsistent with the full deficit of 
halos. 

1.4. This Paper 

The general success of bottom-up galaxy formation sce- 
narios in predicting large-scale structure provides consider- 
able justification for accepting the notion that halos of sub- 
galactic size were already in place at the epoch of reion- 
ization with their baryons somewhat relaxed into their ha- 
los. Yet the paucity of mini-galaxies today suggests that 
cooling could not generally have proceeded to the point of 
star formation, and that somewhat larger clouds, formed 
later but with the greater cooling rates afforded by their 
stronger self-gravity, formed the first stars. 

It is reasonable to suppose that follo wing reionization, 
baryo ns were heated and flowed outward ( Barkana fc Loeb] 
1999), but that they may still be in the proximity of those 
halos. The dark halos of clouds provide a means to restrain 
the dispersal of baryons after the epoch of reionization. 
Though the ability to retain the outflow is dependent on 
the halo model, this could help to explain the apparent 



endurance of low equivalent width clouds, as suggested 
by the low observed v alues of their evolution parameter 
(Weymann et al., 1998). 

High-sensitivity observations of low-redshift Lya clouds 
may provide a means to corroborate these assertions. PSSI 
used HST/GHRS/G160M to study the absorption spectra 
of bright, low redshift Seyfert and BL Lac galaxies at high 
spectral resolution (~ 19 kms -1 , or 80 mA). Their ob- 
servations were supplemented by archival data from the 
same instrument. Some of these spectra reached a "sen- 
sitivity" (see §4.2) as low as about 12 mA. One of the 
important ou tcomes of their analysis (Penton, Shull, & 



Stocke 2000a, hereafter PSSII) was the determination of 



the mean HI equivalent width distribution function /(W). 
Results showed that the redshift frequency of Lya clouds 
with W J£ 12 mA (HI columns N HI £ 10 123 cm -2 for 
b ~ 25 kms -1 ) approach dM/dz ~ 450 per unit red- 
shift. In their modeling of the clouds, PSSII suggested 
that 20-45% of all baryons in the low redshift universe are 
associated with these clouds. 

In contrast to expectations b ased on the simul ations 
mentioned above, both PSSI and Tripp ct al. (1998 ) show 
that significant numbers of clouds (up to ~ 50%) appear 
to be very isolated. One wonders, what exactly is going 
on in galaxy voids? This paper is an attempt to answer 
that question. The paper begins with an analysis of the 
observational basis for various measurable cloud quantities 
(§2), seeks the proper way to specify the environment of 
a cloud, and then shows how galaxy redshift catalogs can 
be used to measure it (§3). The source and nature of the 
data used in this study are presented in §4. Details on 
the method by which mass is attributed to galaxies, and 
how tidal fields are summed, arc presented in §5. Section 
6 presents the distribution functions of clouds specified by 
tidal field upper and lower limits. In addition, a varia- 
tional approach is used to investigate the effects of various 
assumptions on the output distribution functions. There 
is a discussion (§7) and conclusions (§8). 

Though the effect of employing different cosmologies on 
the results of this analysis is minimal, a cosmology must 
be specified, as future modeling will require a cosmological 
context. A flat FRW cosmology with A = 0.7, and h = 
0.75 is used, unless otherwise noted. The cosmic baryon 
density is assumed to be f2&/f2 m = 0.11. 

2. ABOUT CLOUDS 

Before studying these clouds in particular, it may be 
helpful to review their basic observational relationships. 
The cumulative line density dhf jdz of evenly distributed 
clouds of comoving number density n and which produce 
a column density of > Njji (or, alternatively, an EW > 
W) at projected distance r p from the center of a spherical 
cloud is, 

— = ^(1 + ,)»_, (1) 

where dl is the element of proper path length, such that 
the comoving distance coordinate element dr = (1 + z) dl, 
and the column density is formally defined by the integral 
of the HI density over a column of length L through the 
cloud at a projected radius r p , 



N. 



HI 



n H i(r > r p ) dL. 



(2) 



We also have, by way of a definition of 7, 

dz ( dz ) ^ ^ 

Solved for n, Eqs. 1 and 3 yield, 

_ / dAA (1 + z)'<' z dz 
\ dz J 7rr|j dZ 



(3) 



(4) 



For a flat universe with non-zero cosmological constant 
dl R Q 



dz 



(1 + z) A /n m (i + z)3 + n A )' 



(5) 



( Scott, D. and Silk ) where i?o = c/Hq. Assume that the 
radius at which a given column density occurs in an indi- 
vidual cloud is 

r p cx(l + z) e , (6) 

then taking the derivative of Eq. 4 with respect to z we 
find, 

ldn_7-2e-3 d 2 l/dz 2 

n dz (1 + z) dl/ dz 



(7) 



The last term is, 



d 2 l/dz 2 
dl/ dz 



1 + 2 



3/2 



(i + n A /(n (i + z ) 3 )) 



Substitute the above into Eq. 6 to find, 



1 Mu- 
ii dz 



4- 



1 + z 



7 + 



-3/2- 



i + r> A /(fi (i + 2) 3 ) 



2e- 2 



Assume that the comoving number density of some popu- 
lation of clouds is constant so that l/n(dn/dz) = 0. Then, 



7 = 2 + 2e - 



3(1 ± ^) 3 Q- 



2 v /o m (i + z) 3 + r! A 



Let us consider the specific case in which there are asso- 
ciated metal lines. It is probable that the radius at which 
metal lines are detected is fairly constant if they require a 
hot halo to maintain them in ionized form. In that case, 
the cross section for absorption is constant, and e = 0. 
With dn/dz — at z = 0, we find the evolution parameter 
for metal-line clouds is 

Imetal = 1-55, 

for the adopted cosmological parameters (§1), consistent 
with the evolution parameter for metal-line HI clouds ob- 
served byWeymann ct aL (1998). However, when we are 
dealing with primordial Lya clouds that are DM-held, the 
radius at which a given column density is observed may 
change. We decompose e into the effects of a changing 
metagalactic flux and the cloud density profile. Assume 
that 



J(z) = J (l + z) i 
p(r) = p r- s . 



(9) 
(10) 



Equation 10 entails a temporally static matter distribu- 
tion for each cloud. This may be unrealistic; small clouds 



may still be expanding, and larger clouds may be recol- 
lapsing. However, analytical modeling at this level of de- 
tail is clearly inferior to hydrodynamic simulations, so we 
forgo this in favor of the simplification that clouds are 
static. Since optically thin gas has a neutral fraction in- 
versely proportional to Jq(z), and directly proportional to 
n H , then tlhi oc r _2<5 /(l + zp. The column density at a 

cloud projected radius r p is then Nhi oc rp ^ 2<5_1 ^/(l + z)^. 
The impact parameter r p at which the cloud has a neutral 
column Nhi varies as, 



r p oc (l + z)-^ 24 " 1 ). 



(11) 



Apparently we may substitute — £ /(26 — 1) for e. If f ~ 3 
or so, as suggested by the work of [Shul! ct al. (1999 ), and 
S ~ 1.75, a somewhat flattened baryon distribution, then 
e = — 1.2, and at z = 0, Eq. 8 yields, 



Id 



-0.85, 



If we can extend the 7 — EW correlation presented in 
Weymann ct al. (1998 ) down to, say, 12 mA, then the 
above is in agreement. 



But what are the predictions of simulations? As Riedi 



ger et al. (1998 ) point out, if the ionizing flux is approx- 



imately constant, then the column density for clouds ex- 
panding with the Hubble flow is Nhi oc n 2 H L oc (1 + z) 5 . 
Therefore, if Jo oc (1 + zp, where £ ~ 3.0, then Nhi oc 
(1 + z) 2 



To summ arize, the analytical approx imations of Dave 



et al. (1999) and Riediger et al. (1998) are based on the 



presumption of homogeneous clouds of negligible mass 
which expand freely but have sharp edges, whereas the 
view presented in this paper contends that the clouds are 
self-gravitating, and have the gas density gradients char 
acteristic of self-gravitating clouds. As suggested by Dave 



et al. (1999), their cloud model can explain the observed 
low 7 for low EW clouds by citing the effects of a decline 
in flux of the ionizing background with time, but to ac- 
complish this f (Eq. 9) must be of order 4.5, rathe r more 
than suggested by the work of Shull et al. (1999). Our 
preference is that cloud density gradients in conjunction 
with a decline in Jq{z) are r esponsible for the low o bserved 
cloud evolution parameter ( Weymann et al., 1998 ). 



3. THE ENVIRONMENTS OF CLOUDS 

Self-gravitating clouds are delicate structures and can 
be easily disrupted. We are searching for one character- 
istic which by its absence allows clouds to abide, while 
in its presence the cloud is destroyed. Intuitively, isola- 
tion corresponds to the lack of that factor which destroys 
clouds. Hence, our parameter deals with the proximity 
to galaxies and groups of galaxies. Heretofore, the degree 
of isolation of clouds has been generally characterized by 
the distance (in projected radius and velocity) to the near- 
est galaxy. However, one megaparsec from a small galaxy 
is not the same as one megaparsec from a giant. Opti- 
mally, one should weight the influences of galaxies by their 
mass and distance. Besides the effects of ionizing radia- 
tion, which is here assumed to be uniform over the small 
range of redshifts of the data to be studied, the principal 
constraints on the physical integrity of clouds is from shear 
forces or ram pressure, and the strength of the tidal field. 
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In addition, the heightened probability that a small halo 
would be accreted to a larger one is a direct function of 
ambient density; it is expected to produce an effect similar 
to that of a tidal held - a flattening of the cloud EW dis- 
tributions in high density areas. For isolated clouds, only 
tidal holds are relevant, though the flux of ionizing photons 
may cause heating which will help to disperse the baryonic 
content of small clouds. Because a uniform distribution of 
mass will have no tidal effect on a cloud, galaxies are the 
main sources of tidal fields on large scales. 

3.1. Tides produced by galaxies 

If one takes the second spatial derivative of the summed 
gravitational potentials of galaxies at a particular point, 
one gets a quantity indistinguishable from the tidal field. 
One can imagine, then, that minima in the magnitude of 
tidal fields would occur in galaxy voids. The tidal field 
T is in fact a function expressible in units of inverse time 
squared, so that if it is multiplied by the Hubble time 
squared it becomes a dimensionless quantity T = T/Hq. 
The dimensionless form T is hereafter referred to as the 
"tidal held" parameter. Tides may be of vector or scalar 
form. In this paper, a scalar tide is calculated, though 
we must first consider its vector form. For an individual 
galaxy, the only non-zero component of T is radial; 



J_ _d_ ( GM{R) \ 
R H 2 dR \ R 2 J' 



(12) 



In summing the effects of many galaxies at a given point 
along the line of sight (LOS), we must decompose this into 
Cartesian components. For instance, the x-component, 
transverse to the LOS, is 



2GM^ 
H%R 4 1 '' 



(13) 



where R — yj Ax 2 + Ay 2 + Az 2 . The components trans- 
verse to the observer's LOS are calculated using the angu- 
lar diameter distance pertinent to the cosmological model 
employed in this paper. The proper displacement Az is 
found using Eq. 5 (§2). That is, the radial (LOS) posi- 
tion is determined by the redshift alone; no allowance is 
made for peculiar velocities. This causes some positional 
error, which will be discussed in §5.5. In this manner, the 
contributions of each galaxy can be added, then combined 
in quadrature to arrive at an averaged scalar tidal field T 
acting on a cloud, 



T = 



T 2 + T 2 + T 2 



(14) 



Let us see how the tidal field works to disrupt clouds. 
The condition for stability to tidal disruption is as follows, 



H 2 Tr < 



Gm(r) 



(15) 



where m(r) is cloud mass summed to radius r. When a 
cloud is close to a single galaxy, then the radial form will 
be more pertinent. By application of Eq. 12 to single 
galaxy, the tidal held is, 



_ GM(R) ( _ d log A/ 
R H 2 R3 \ d logR 



(16) 



For an isothermal halo, d log M/d logi? is unity when 
R < R t , a nd zero beyond that. The standard h alo model 



is that of |Navarro, Frenk, fc White 1996j , |1997| (hereafter 
NFW), which we discuss in more detail below (§5.1.2). 
The above derivative for the NFW halo does not approx- 
imate a constant function of r unless baryons "conspire" 
to produce a flat rotation curve inside r max . However, be- 
yond r max the density profile tends toward an index of —3, 
and a simple solution for Tr no longer maintains. 

The simplest halo is isothermal, possibly truncated at 
some characteristic radius Rt- The mass of such a halo 
increases linearly with radius, so that one may describe 
this radial mass distribution by a constant, JC = M(R)/ R. 
The circular velocity of such a system is v c = \J G7C, and 
the total mass is M to t — ICRt ■ F° r a truncated isothermal 
halo, Eq. 13 yields, 



1 R — \ l&Mtot 

imr 



.(R > R t ). 



(17) 



Note that the effect of halo mass truncation on the tidal 
field is to produce a discontinuity at Rt- While a sharp 
truncation is not realistic, it will not produce a significant 
effect on the distributions of clouds as we shall in all cases 
assume a total mass. 

In the context of a specific cloud, Eq. 15 can be written, 



G 



< 



4tt 



p(r) 



(18) 



where p(r) is the average over-density of the cloud within 
radius r. Thus, the requirement for stability can be stated, 



P> 



ZH 2 T 
AttG 



2T Pc 



(19) 



where p cr it is the critical density, p cr u — 3Hq /8irG. There 
is one important caveat to this: the proper number den- 
sity of galaxies declines with decreasing redshift due to the 
expansion of the universe. Therefore, since the average dis- 
tance of galaxies from an average void cloud is increasing 
with time, it is clear that in the past the tidal field on 
a void cloud varies roughly as T cx (1 + z) 3 . However, 
clouds which are in some proximity to a galaxy may have 
a weaker T with lookback time if the peculiar velocity to- 
ward the galaxy is greater than the HR. When we think 
of the tidal field affecting void clouds we must antedate 
to some extent the current T to get a feel for what it has 
been subjected to in the past. For instance, at a lookback 
time of 7.2 Gyr [z = 1), the tidal field affecting a void 
cloud will be ~ 8 times the present value. 

Self-gravitating clouds have negative radial density gra- 
dients. The p > 27" Pcrit criterion applies differentially to 
the clouds within which M(r)/V(r) > 2T p cr n- Though 
tidal fields affect all clouds, only the smallest are de- 
stroyed. Thus the effects of tides will be manifested over 
the distribution of cloud sizes, and hence over the whole 
range of column densities and EWs. The physical effect 
will be a tidal truncation, reducing the average column 
density a cloud would produce for a given halo velocity. 
It is expected to affect low-mass clouds more dramatically 
than more massive clouds. 

A cloud near a galaxy will experience a radial tidal field, 
and it is unlikely to have experienced a higher tide in the 
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past, given its probable infall into the galaxy's (or group) 
potential well. An example may clarify the utility of Tr. 
A homogeneous cloud with an average over-density p = 
%Pcrit would be marginally destabilized by a tidal field of 
Tr = 1.0. A galaxy with a flat rotation curve of v c = 220 
km/s would destabilize this cloud at a distance of 2.08 
Mpc if the galaxy halo truncates at 0.5 Mpc, or at 2.94 
hj§ Mpc if it has an untruncated isothermal halo. Thus 
we would not expect to find discrete clouds of low EW in 
close proximity (<; 2 Mpc) to galaxies. A galaxy with an 
NFW halo of the same velocity would have a tidal field 
T = 1.0 at a distance slightly less than 2/3 times that of 
the truncated isothermal halo. 

We now turn our attention to the problem of sorting 
clouds by their environments. If the tidal field is calcu- 
lated along the lines of sight to AGN, then if one specifies 
a limiting tide T, absorption systems can be separated 
into two catalogs, one with T providing the upper limit, 
denoted C(T U ), and the other with T as a lower bound, 
C(T L ). The actual calculated value of the tidal field will 
depend on which galaxy halo model is used. The NFW 
halo is the outcome of analysis of N-body simulations, and 
has a "cuspy" core approaching a density p cx r _1 . How- 
ever, more recent, higher-resolution simulations find the 
core of the halo t o be even more cus py than the NFW 
model (p cx r -15 ) (Moore et al., 1999b), suggesting a con- 
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Burk 



flict with observed galaxie s (Firmani ct al 
ert, 2@; |Wu k Xue, 2000| ). Furthermore, the NFW halo 
may be inconsistent with what is suggested by the steadi- 
ness of the radial dependency of the velocity dispersions 
of dwarf galaxies about their parent giant field galaxies to 



500 h^l kpc ( |Zaritsky fc White, 1994] ; |Zaritsky 



radii r. 

et al., 19971), far beyond the radius at which the halo has 
a maximum circular velocity (r max is of order 15-40 kpc). 
Because of these perceived problems with the NFW model, 
and the simplicity of the truncated isothermal halo, the 
latter is used as a standard model in this paper, relegating 
the NFW to the position of first alternative model. 

It is possible to coarsely see the effect of tides in the 
original PSSI data by binning absorption systems accord- 
ing to the projected radii (with Av < 500 kms" 1 ) to 
the nearest galaxy. The results show that very low col- 
umn density systems are not seen within ~ 1 Mpc of a 
galaxy (Fig. 1, top panel), but the farther one goes from 
the nearest galaxy (middle and bottom panels), the lower 
the Doppler parameters and column densities which are 
recorded. If void clouds were freely expanding, one would 
not see lower b- values with decreasing EWs, as the bulk 
velocity of expansion is added to the thermal broadening. 

These provisional findings are consistent with the no- 
tion that tidal fields disrupt smaller clouds, and also sup- 
ports the concept that the expansion of baryonic clouds 
are at least restrained by the gravity of their dark halos. 
It is reasonable to suppose that outside of massive halos, 
tidal fields provide the main constraint on cloud stability. 
We therefore formally adopt the tidal field T as the dis- 
criminator for cloud environments; a low T then indicates 
an environment in which cloud density is essentially un- 
constrained, and a high T would mean that clouds with 
density p r j <J 2Tp cr i t would be disrupted. 

4. THE DATA 
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Fig. 1. — Column density as a function of b- value for data from 
PSSI for 3 bins of projected radius to the nearest galaxy. Radius was 
derived using a "retarded" Hubble flow model (clouds are placed 
at the galaxy redshift when they are within ±300 kms -1 of the 
galaxy). Cloud b- values appear to have lower limits which correlate 
inversely with galactocentric radius. The variation of the lower lim- 
its of column densities with radius are consistent with the effects of 
tidal fields. 



The basic data (PSSI) required for calculating the distri- 
bution /(W) as a function of tidal field are of three types; 
cloud absorption-line data (wavelength A, rest equiva- 
lent width W, Doppler parameter b), sensitivity func- 
tions, and galaxy redshift catalogs. We begin with the 
HST/GHRS/G160M data of PSSI along 15 sightlines to- 
ward luminous, low-redshift Seyfert and BL Lac galaxies. 
Table one presents 15 target RA and Dec as well as the 
spectral range covered by the data for each line of sight. 

Table 1 

Spectral targets and wavelength range 



Target 


RA 


Dec (.T2000) 


Xmin (A) 


X„«« (A) 


3CJ 273 


12:29:06.7 


+02:03:08.6 


1214.11 


1301.46 


Arkclian 120 


05:16:11.4 


-00:08:59.0 


1222.51 


1258.75 


ESO 141-G55 


19:21:14.3 


-58:40:13.0 


1222.51 


1265.92 


Fairall 9 


01:23:45.8 


-58:48:20.0 


1219.79 


1276.38 


H1821+643 


18:21:57.3 


+64:20:36.4 


1231.66 


1267.75 


IZW 1 


00:53:34.9 


+ 12:41:36.0 


1221.40 


1257.54 


Mrk 279 


13:53:03.4 


+69:18:29.6 


1222.53 


1258.74 


Mrk 290 


15:35:52.4 


+57:54:09.2 


1231.76 


1268.99 


Mrk 335 


00:06:19.5 


+20:12:10.5 


1221.41 


1257.55 


Mrk 421 


11:04:27.3 


+38:12:31.8 


1221.43 


1257.57 


Mrk 501 


16:53:52.2 


+39:45:36.6 


1221.42 


1257.56 


Mrk 509 


20:44:09.7 


-10:43:25.0 


1219.47 


1267.84 


Mrk 817 


14:36:22.1 


+58:47:39.4 


1222.57 


1258.77 


PKS 2155-304 


21:58:52.1 


-30:13:32.1 


1222.58 


1293.67 


Q1230+0155 


12:30:49.9 


+01:15:23.0 


1216.96 


1254.23 



4.1. The Absorption Systems 

Source spectra probe sightlines to low-redshift active 
galaxies brighter than V < 14.5 mag (Table 1). Only two 
of the sightlines include data beyond a redshift z = 0.045. 
The well-studied sightline to 3C 273 includes clouds with 
the largest maximum redshift of z ~ 0.07. 
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The data include cloud wavelength and error, equiva- 
lent width, and its error, the resolution-corrected Dopplcr 
parameter, and its error. Some of the spectra are pre-, 
and some post-deployment of COSTAR (Corrective Op- 
tics Space Telescope Axial Replacement), but there is no 
attempt to separate them except insofar as their sensi- 
tivity functions are calculated somewhat differently. No 
clear trend in Doppler parameter is noticed between pre- 
and post-COSTAR data. 

The clouds are considered to be randomly selected from 
a parent population - one which is effected by ambient 
tidal fields - within the constraints of the sensitivity of 
the instrument along the line of sight (LOS). The reader 
is referred to PSSI for a detailed description of the reduc- 
tion of the HST data. Their procedure involves fitting 
the continuum to a polynomial by increasing the order of 
the polynomial until no further reductions in x 2 are possi- 
ble. All negative spectral fluctuations from the continuum 
greater than 1-tr are treated as the original absorption line 
list, and the lines are fitted with gaussian components ac- 
cepting only those with values 12 < b b s < 100 kms -1 . 
The significance levels of clouds are calculated by integrat- 
ing the signal to noise ratio per resolution element over the 
cloud (see below). 

4.2. Sensitivity functions 

The sensitivity functions indicate what cloud EWs can 
be observed given the significance level imposed for the 
cloud catalog. 

Let us begin with a review of the relationships between 
signal-to- noise (S/N), significance level (SL), EW, and 
sensitivity (<S). The signal S of an absorption system is 
the integral of the counts-deficit relative to the continuum 
over a resolution element, 



(fcontW-fx)S G (t/4)d\ = An, 



(20) 



RE 



where Sq is the GHRS sensitivity at A in units 
counts s _1 diode -1 ergs -1 s -1 cm -2 A -1 (GHRS handbook), 
and the absorption system is represented by the difference 
of the continuum f con t and the observed flux f\. The 
time t is divided by 4 because the exposures are "quar- 
ter stepped" using FP-SPLIT to help remove the effects 
of photocathode granularity and diode to diode detector 
variations. The noise N , based on the total counts in the 
continuum integrated over one RE, is 



N 



Ire 



fcont S G {t/A)d\ = s/n, 



so within a given RE, 



S/N=^. 
\/n 



(21) 



(22) 



Other sources of noise are said to be small so that the S/N 
is within <~ 5% of the "root n" photon statistics (Final 
Report, GHRS Science Verification Program, 1992). It will 
be assumed that "root n" statistics are accurate though 
it may be a small under-estimate of the true sensitivity 
in mA. The resolution element RE is estimated at 1.2 
diodes for pre-COSTAR and 1.1 diodes for post-COSTAR 



installation (PSSI). The observed equivalent width of a 
cloud within a resolution element is, 



EW 



An 
n/RE 



An 



-RE, 



(23) 



where RE ~ 80 mA(PSSI). 

A cloud detected at a 4.0-ct SL in one RE has an EW 
W = An RE I n such that An = 4.0\/n- The sensitivity 
per RE, S, is the EW of a 4-cr cloud. Thus, 



n y/n 



(24) 



These sensitivity functions are plotted as solid lines in Fig. 
2. 

The total cloud EW is the sum of the EWs of each reso- 
lution clement. The total significance level of the cloud is 
the imposed significance level times the summed equiva- 
lent width per RE divided by the sensitivity per RE. Using 
the equations above, terms cancel, and the SL of the cloud 
is, 



SL = J2 



An(i) 



(25) 



Thus the SL of the cloud is the sum of the S/N in each 
RE of the cloud. In practice this is done by integrating the 
S/N over the best fit Gaussian. Thus, if a cloud of a given 
EW is found in one RE, and if its EW > S, then a cloud 
with a similar EW, but spread over more RE's, will have 
the same significance level, for the counts deficit per RE 
An oc b~ x , while the number of RE's spanned by the cloud 
is oc b; in combining them the b— factors cancel. Thus, if 
a cloud appears above the sensitivity function, then it is a 
4-sigma cloud, regardless of its b- value. 



..(a) 



1 1 1 r 

3C 273 



H 1 1 1 1 1 H 



(b) 



PKS 2155 t / 



1220 1240 1260 1280 1300 

Wavelength A(A) 

Fig. 2. — The sensitivity functions (panels a-o) as derived by the 
author from the flux (solid line) and extracted from PSSI on the 
basis of fluctuations around the continuum. The dashed lines are 
the sensitivity functions from PSSI. Clouds with significance greater 
than 4-cr (as assigned in PSSI) are represented by solid dots; open 
dots are for clouds with SL in the range, 3-cr < SL <4-<j 

The sensitivity functions published by PSSI (the dashed 
lines in Fig. 2) are based on the noise in the spectrum 
around the placed continuum. Clouds are plotted at their 
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Fig. 2.— (Continued) 
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Fig. 2.— (Continued) 



observed Lya wavelength, and rest EW. Filled circles rep- 
resent the PSSI assignment of SL values > 4 — a, while 
clouds with 3-er < SL < 4-c are shown by open symbols. 

It is apparent that the two sets of sensitivity func- 
tions do not always agree. While many are in substantial 
agreement (e.g., Mrk 421, Mrk 279, Akn 120, Mrk 817, 
Q1230+0115, ESO 141-G55 and possibly Mrk 509), seri- 
ous problems are apparent for IZW 1, Mrk 290, Mrk 501, 
Fairall 9, and H1821+643. Some sightlines appear to have 
clouds with SL > 4 — a and EWs less than the sensitivity 
function, or clouds with SL < 4 — a with EWs larger than 
the sensitivity function (e.g., Figs. 2a, d, e, j, k, and m), 
which violates the derivation given above for a 4-cr sensi- 
tivity function. This may indicate a source of error which 
is outside the normal physical explanations. There may be 
"grey" areas in the data processing which can introduce 



substantial errors. For instance, unrecognized 3-cr clouds, 
and possible continuum fitting errors may cause an over- 
estimation of the continuum noise, and a correspondingly 
under-estimated significance level of a given cloud. There 
are four cases in which cloud EWs are significantly larger 
(in mA) than the flux-based sensitivity function at that 
location (3C 273, IZW1, Mrk 335, and Mrk 501) which 
may be explained by this logic. But there are a similar 
number in which their sensitivity functions are lower, and 
this cannot be explained in the same manner. The cause of 
this variation is uncertain. However, since the sensitivity 
functions based on the flux in the continuum is a simpler 
relation, the author judges that the flux-based sensitivity 
functions are more likely to reliably reflect the actual sen- 
sitivity. For this reason our standard processing of data 
uses flux-based sensitivity functions. The significance lev- 
els of the cloud are hence re-evaluated according to their 
EW relative to the sensitivity function. 

4.3. The Galaxy Catalogs 

A catalogue of galaxies is constructed along each LOS, 
and this is used this to calculate T. All galaxy catalogs 
are taken f rom the CfA Redshift catalog, maintained by 
J. Huchra (Huchra ct al. 1990], [Huchr a, Geller, fc Corwiri 



1995|, IMarzke Huchra, fc Geller 1996|, Progin Gellcr, & 



Huchra 1998| , |Huchra, Vogclcy, fc Geller 1999| )|j. The data 
which the catalog has that is of use to our enterprise is the 
position (B1950), velocity (with error), magnitude (with 
error), and the de Vaucouleurs T-Typc. The T-Type is 
negative for early-type galaxies and quasars and positive 
for spirals and irregulars. This allows one to pick out those 
galaxies which are not well-described by the Tully-Fishcr 
relation. The catalogs specific to each line of sight include 
all galaxies in the CfA catalog which are within a cylinder 
7.5 h^ 1 Mpc in radius, and within 15 degrees of the line 
of sight, and which have a redshift range extending 750 
kms^ 1 beyond the spectral endpoints (see Table 1). We 
discuss the completeness levels of these catalogs in §5.4. 

The radial positions of both galaxies and clouds relative 
to the observer are assumed to be that attributed on the 
basis of the redshift alone; clouds and galaxies are assumed 
to have a vanishing peculiar velocity (see §5.5 for further 
comment on this). Galaxy apparent magnitudes are as- 
sumed to be blue magnitudes, though the specific filter 
varies, and there is a low level of "contamination" with 
other filters when blue is not available. The Tully-Fisher 



relation (Tully & Fisher 1977 , hearafter, TF) is used to de- 
termine the circular velocity range (see §5.1.1), and from 
this, the galaxy mass as a function of halo model. Though 
not all galaxies are spirals, we can adjust for this by noting 
that galaxies with negative de Vaucouleurs T- Types gen- 
erally have mass-to-light ratios 2 to 4 times greater than 
spirals. We discuss this detail in §.6.2. 

Much depends on the reliability of the apparent mag- 
nitudes in the assembled galaxy catalogs. For those in- 
terested in this level of detail, the Appendix offers a rec- 
onciliation of the CfA magnitudes with magnitudes ex- 
tracted electronically from digitized Palomar Observatory 
Sky Survey plates, where we conclude that the CfA mag- 
nitudes are accurate to ±0.4 mag.. 



see http:/ /cfa-www. harvard. edu~huchra/ 
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5. THE GALAXY MASS AND TIDAL FIELDS 

The next task is to calculate the tidal field produced by 
the galaxies in each catalog for each point along the lines 
of sight. We must first discuss the derivation of galaxy 
masses, then the tidal field. We also inspect the results 
of applying this methodology to the galaxy catalogs along 
the sightlines to AGN. 

5.1. Galaxy mass 

The basic procedure for deriving a galaxy's mass is as 
follows: apparent magnitude and redshift are used to de- 
rive an absolute magnitude, then an empirical relation 
must bo exploited to ootimato the galaxy maaa. The over — 
limi 



whel ming fraction of galaxies are of spiral type . For those, 



the TF relation is used to determine the circular velocity of 
the galaxy, and hence its mass as a function of the adopted 
halo model. The TF relation expresses the luminosity as 
a power of the circular velocity, 



(26) 



where v c is the maximum circular velocity of th e galaxy, 
and the slop e (3 is a function of spectral range. Tully & 
Pierce (2000| ) show that the blue-band slope, relevant to 
the CfA catalogue, is (3b — 2.91. The B-band galaxy lu- 
minosity is influenced by recent star formation, and so is 
more likely to experience larger excursions from the mean 
relation than, for instance, the K-band relation. For in- 
dividual galaxies, there is a 0.55 mag dispersion relativ e 
to the calibrated B— band relation (Tully & Pierce, "2000| ), 
but only 0.4 mag in R and K' . Given a DM halo model, 
one can extract a galaxy mass with error. 

The derived mass is a function of the halo model. But 
having accurate relative masses of galaxies is more impor- 
tant at this stage than achieving accurate absolute masses. 
We consider two halo types. Our standard model is a trun- 
cated isothermal halo, and the NFW model, is treated as 
the first alternative. We treat the two halo models in order 
below. 

5.1.1. The mass of truncated isothermal spheres 

We first assume that the distribution of galaxy mass is 
well-modeled by a truncated isothermal sphere. In this 
case, the halo mass density can be approximated by, 



P {R) 



K 



AttR 2 1 ' 



■(R < Rt), 



(27) 



where K. is called the total mass distribution constant, and 
R t is the mass-truncation radius. Since the TF relation is 
based on circular velocities for dark matter plus baryons, 
K, is taken to be the observed mass-distribution constant, 
where 



K 



G 



(28) 



Assume that the galaxy masses are truncated at a radius 
such that the density at R t is the same for all galaxies. 
A glance at Fq. 97 will show that the truncation radius 



mce a 

Rt oc 1C 1/2 . With the above relations (Eqs 26, 27, and 28) 
we may construct the following scaling relations valid for 
isothermal halos: 



v* 



L_ 

C* 



V0f 



K_ 

~K* 



R4_ 

R* t 



( M{R t 



\M*{Rt) 



(29) 



This halo model is standardized by requiring that an C* 
galaxy has a truncation radius Rt of 0.5 h^ 1 Mpc, con- 



sist ent with the observat ions by Zaritsky & White (1994) 
and Zaritsky et al. (1997) of the steadiness of velocity dis- 
persions of satellite galaxies at large galactocentric radii 
around parent field spirals. 

With the use of the cataloged redshift, the absolute mag- 
nitude of the galaxies is calculated. The Tully-Fisher re- 
lation is used to derive the circular velocity. The relevant 
error in absolute magnitude is the CfA magnitude error 
(±0.4 mag, Appendix) combined in quadrature with the 
TF error per galaxy (0.55 mag in the £?-band), which re- 
sults in a tot al error of 0.68 mag. According to Tully & 



Pierce (2000), the maximum circular velocity of a galaxy 



with absolute magnitude of Mb is, 



158 exp 



-(M B + 20.11 -5 log /t 75 ) 
7.27 



+0.255i> c 



-0.19-Uc 



kms 

where the l-er errors in attrib uted total cir cular velocities 
are shown. The calibration of Tully (1988]) indicates that 
M* B = -20.18 
galaxy is then, 



M B = -20.18 + 51ogh 75 . The total mass M* of an C* 



AT 



G 



(31) 



where v* is the circular velocity of a C* galaxy. With the 
scaling equations (Eq. 29) we see that M oc vf. The mass 
of a galaxy may be expressed as, 



M = 



M I — 

C* 



-i +0.54Af 



J -0.35M 



(32) 



where L/C* = C exp(-0.4M B ), and C = 8.47 x 10" 9 , 
and (3b = 2.91. The 1-a range in observed mass M is then 
0.65M to 1.54M. 

Not e that if one sets M B = M% = -20.18 mag (h = 
0.75) ( fTully, 1988D , one can solve for v* by substituting 
M B and using Eq. 30. In the B-band, 



161.5 kms 



(33) 



What evidence, you might ask, is there to support this 
halo model (other than flat rotation curves)? Each halo 
at large galactocentric radius, and the gas associated with 
the halo, is assumed to scale as shown in Eq. 29, indicating 
that at a radius which scales r oc C 1 ^ , the gas density of 
each system would be the same. We assume that the EW 
of HI gas found at distance r p from a galaxy should be 
proportional to th e gas density there. A recent paper by 
Chen et al. (2001) has calculated the radius at which a 
characteristic absorption of 0.3 A occurs and found that 
in the K-band, r oc £^ 28±0 08 a t a characteristic radius 
r ch fa 180 h- 1 kpc. In the B-band, r oc £^ 39±0 09 . The 
isothermal halo model predicts that the exponent for both 
is the respec tive which, given that (3k =3.51 (Cole 
et al, 200 jj), and (3 B = 2.91 yields r oc £ 285 , and oc 



£ U J44 , both well within the er rors quoted above . That is, 
the gaseous halos observed by Chen et al. (2001), and the 
scaling laws for isothermal halos, are consistent with each 
other. 

On the other hand, the truncated isothermal halos pro- 
posed here are significantly more massive than those of 
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the favored, NFW halo. While the mass within the dis- 
tances probed in studies of the ma ss-luminosity ratio (e.g.. 
Bahcall, Lubin, fc Dorman, 1995 ) is compatable with the 
NFW halo, so is the mass of an isothermal halo. It is in 
the realm beyond the virial radius, which has not been 
effectively probed, where this extra mass may be found. 
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Fig. 3. — The ratio of the NFW virial mass to an isothermal 
halo truncated at a radius scaled to that of an C galaxy with n = 
0.5 Mpc. The assumption of an isothermal halo will have the 
effect of attributing relatively more mass to small halos than to 
large. The skew from strict proportionality will be less than about 
10% over the range of halo velocities of interest in this paper. 



5.1.2. The mass of a NFW halo 

If halos are described as NFW, then /C = M max /R max , 
where R m ax is the distance at which the circular velocity 
of the halo is a maximum. The density is, 



p{R) = 



Pcrit 



(R/R s )(l + R/Rs 



(34) 



(NFW, 1997), where R s is the scaling radius; R ma x — 
2.16i? s . NFW halos are defined for the dark component 
only, and one has to adjust for the baryon fraction. This 
can be done by requiring that enough baryons are added 
within r max to give a flat rotation curve inward of that 
point. As with the isothermal halo, the circular velocity 
is derived from the apparent magnitude and the redshift 
(refer to Eq. 30). The concentration parameter c is defined 
as the ratio of the virial radius -R200 and the scaling radius 
Rs- 

By fitting the relation in Fig. 9 of NFW (1996) we find, 



«200 



= 2.44u 



-0.110 



Figure 10 in NFW (1996) gives, 

Ft - 1 5 x in -2 ?; 1 - 49 



(35) 



(36) 



From Eq. 3 of NFW (1997) we find a reasonable fit yield- 
ing, 



.1/2 



^200 (2.16{^(l + c)-c/(l + c)}) 1 /2 



1.27 c 



,0.28 



(37) 



Equating Eqs. 39 and 41 yields 
c= 132.3 v. 

Next we find, 
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Using these equations we can construct scaling relations 
analogous to those of Eq. 29, 
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The errors propagate in the same way as with the isother- 
mal halo though the virial mass is proportional to oj,^, 
rather than v^ ax for the isothermal halo. This accounts 
for the rising trend in the ratio of masses of NFW and 
isothermal halos of the same v max (see Fig. 3). 

How does this halo accord wi th the characteris tic ra- 
dius for absorption observed by Chen et ah (2001)? 



We 

would expect that it should scale with the virial radius. 
From the scaling relations above, i?200 oc £ 1 /(°- 76 ' 3 ) j yield- 



ing R 200 (K) cx C° K 375 , and R 2 oo(B) cx ££ 4&2 , which are 
supposed to compare with ex ponents 0.28 ± 0.0 8, and 
0.39 ±0.09, respectively from phen et al. (200l|). The 



predicted .PC-band relation is outside the error bars, while 
the £?-band relation is within the 1-cr error bars, but near 
the edge. If instead one were to require the characteristic 
radius to scale with m ax, then the lack of concordance is 
exascerbated, giving exponents 0.43 and 0.52 for K and 
i?-bands, respecti vely, both well-outside 1-cr error bars of 
Chen et al. (200l|). 



The NFW halo, however, is defined by the dark compo- 
nent only. By adding baryons until a flat rotation curve is 
maintained out to r max , one finds that the circular veloc- 
ity is increased by from ~ 12%, for small clouds, to ~ 8% 
for large galaxies. The NFW halos are corrected for their 
baryonic part when comparing the masses of NFW and 
isothermal halos. 

5.2. Error in Attributed Tidal Fields 

We have seen in general terms how the tidal field is cal- 
culated along the LOS (§3.1). If we assign no error for 
the distance R(i), the fractional error in individual tidal 
fields is the same as that in the attribution of mass to a 
galaxy on the basis of its redshift and apparent magnitude 
raj - that is, on order 50%. The error in R(i) will al- 
most entirely result from positional errors along the line 
of sight (see §3.1). The fractional error in the attributed 
tide produced by a number of different galaxies will gener- 
ally be significantly less than from a single galaxy because 
the individual errors are independent. However, since the 
masses and tides vary, there is no simple way to calculate 
the total error except by a specific convolving with redshift 
and mass errors of the effects of each galaxy. This level 
of sophistication is not pursued here, however, we will dis- 
cuss the errors resulting from cloud peculiar velocities and 
redshift errors (§5.5). 

Another source of error is the volume from which the 
galaxy catalogs are drawn. Let us consider the tidal field 
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of an C* galaxy at the edge of the 7.5 Mpc cylinder at the 
same redshift of a cloud (a minimum maximum distance, 
so to speak). For an NFW halo the vector tide at 7.5 
Mpc is Tr ~ 0.00195, while the isothermal halo is ~ 5 
times greater, Tr ~ 0.0098. This tide would produce a 
scalar tide of T = 6.5 x 10~ 4 , and 3.3 x 10" 3 for the NFW 
and isothermal halos, respectively. Thus our estimated 
relative tidal fields may begin to produce serious errors at 
respective tides of this order. We shall see (§5.4, Fig. 5) 
that only a small fraction of galaxy mass is contained in 
galaxies brighter than L* , so the magnitude of the error is 
generally small. 




Fig. 4. — Tidal fields along the lines of sight to target AGNs. Pan- 
els a and b show 3c273 and PKS 2155-304. The tides are expressed 
in units of the tidal field parameter T (§3.1), which is convenient 
for constraining average cloud density for stability against disrup- 
tion (see Eq. 19). The solid lines represent all galaxy data, while the 
dotted lines represent the effects of only those galaxies brighter than 
Mb = —19.8 mag, the absolute magnitude of a tub = 16 galaxy 
at z = 0.036 (see §5.2). The highly significant (4-cr) absorption 
systems, in terms of the author's sensitivity functions (solid lines, 
Fig. 2 a-o, are denoted by the solid circle), while the 3-cr clouds are 
represented by open symbols. The cloud EWs are denoted by the 
error bars, with length proportional to logW. The vertical dotted 
lines show the limits of the useful portion of the absorption system 
surveys (see §3.1 and Table 1). 
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Fig. 4. — (Continued) 
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Fig. 4. — Continued) 
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5.3. Sightline Characteristics 

As mentioned above (§5.2), the LOS component of the 
tidal field is calculated assuming a pure Hubble flow of 
galaxies, while those normal to the LOS are calculated 
using the appropriate angular diameter distance relation. 
Figure 4 shows the calculated tidal fields (solid line) as 
a function of position along the 15 LOS. Cloud positions 
are noted by the circles (filled, SL > 4-cr, open, SL < 4- 
&) with vertical bars, whose lengths are proportional to 
log W. Note the similarity of the tidal fields in the sight- 
lines to 3C 273 and Q1230+0155, which both pass through 
the Virgo cluster. Some stretches of space in the various 
lines of sight have no cataloged galaxies within the search 
criteria, and so have a formally zero tide. A value of 0.002 
is added to all tides so that the line would appear on the 
plots. The significance of the dotted lines which "shadow" 
the tidal field traces is discussed in the next subsection. 



5.4. Effect of Lack of Completeness 

Though galaxy detections are surface brightness se- 
lected, the decision to extract a redshift and place an ob- 
ject in a galaxy catalog are often determined by the total 
apparent magnitude. The CfA catalog is generally based 
on Zwicky magnitudes and have a completeness limit of or- 
der 15.5 < m < 16.0 mag, but the catalog accepts galaxy 
data from a wide range of sources. To get an idea of how 
much mass is unrepresented in the catalogs, a fractional 
cumulative mass function is calculated. Figure 5 shows 
the cumulative fractional mass function for the truncated 
isothermal halo in terms of the absolute magnitude of the 
galaxy with a faint-end slope a = —1.2 (solid line). It is 
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produced by the following integral, 

JJ C M(v c )4>{v c )dv c 



/m(< v c ) 



C f M{v c )^{v c )dv c 



(41) 



where Vi = 37 kms _1 and Vf — 375 kms" 1 , which encom- 
passes the overwhelming majority of galaxian mass, and 
4>(v) is the B-band luminosity function transformed to a 
velocity function by use of the TF relation. With the use 
of Eq. 30 we express halo velocity in terms of absolute 
magnitude. Note that ~ 35% of the galaxy mass is in the 
r ange Mb < —19-8 mag. A galaxy with Mb = —19.8 is 



0.74 L '* (v n ~ 146 kms ), and has an apparent magni- 



tude of 16 mag at z ~ 0.036 

However, it may be reasonable to suggest that the sub- 
halos within a given halo are characteristic of the system, 
and not an excess over the halo attributed to the luminos- 
ity of the central galaxy. The NFW halo is the p rofile for 



an in dividual halo . Studies of "secondary" i nfall ( Zaroub 
& Ho fman, 1993 ; Avila-Rcese et al., 1999| ; Klypin ct al 



1999^ suggest that sub-halos are not included in the mass 
distribution of the standard NFW halo. The truncated 
isothermal halo, however, describes a system extending far 
beyond the galaxy which has formed at its center. In this 
case, the sub-halos are interpreted as an aspect of the sys- 
tem referred to by the mass distribution constant K, (see 
Eq. 28). So in the case of the isothermal halo, we count 
sub-halos twice unless a factor representing the fraction of 
halos with a given halo velocity which are not sub-halos 
is inserted into both upper and lower integrals of Eq. 41. 
The fraction of ha los of velocity v r which a re isolated has 
been estimated by Avila-Reese et al. (1999 ) , and at low v c 
it appears to follow the relation /j so = 1.36(w c /225) — 0.23 
for halos smaller than 158 kms -1 (0.7 x 225). The re- 
sults of applying /j so to the the numerator and denom- 
inator of Eq. 41 is the dotted line in Fig. 5. While 
with the former Mb = —19.8 mag occurs at Jm — 0.65, 
with the latter it occurs at fm = 0.44. The implication 
of this is that the sub-halos of the brighter galaxies are 
counted already, so less is missed. On the other hand, 
it now appears possible that we may have over-estimated 
the tidal field at very low-redshift by observing their sub- 
halos, and thereby effectively counting them twice. There 
is thus an inevitable tendency to find more mass at low 
redshift. Above z ~ 0.036, about half of the mass in galax- 
ies is missed for catalogs with a uniform magnitude limit 

TUB = 16. 

Four spectra extend above z = 0.036, 3C 273, PKS 2155- 
304, Fairall 9, and H1821+643 with a total of 21 clouds. 
To test for the effects of incompleteness one could calculate 
the equivalent width distribution with, and without, the 
data from z > 0.036. If there is no significant difference, 
then one might say that completeness issues are not debili- 
tory to the conclusions of the paper, for if no completeness 
effects found by the removal of data with z > 0.036 then it 
is unlikely to do so at lower redshifts. However, we should 
expect that clouds at higher redshift may have calculated 
ambient tidal fields a factor of ~ 2 low, hence one runs the 
risk of placing what may be GDS clouds in void catalogs. 

A sense of the potential effects of incompleteness can 
be seen by placing a uniform limit on the luminosity of 
galaxies used to calculate tides. In Fig. 4 (panels a-o) the 
dotted line represents the tidal field when galaxies with 



Mb fainter than —19.8 mag are removed from considera- 
tion. Generally the discrepancy between the total and the 
contribution from only the brighter galaxies declines with 
increased redshift as fewer and fewer low luminosity galax- 
ies make the catalog (a good example of this is panel b). 
This is characteristic of fields of roughly uniform apparent 
magnitude limits. On the other hand, the sightline to 3C 
273 (panel a) has regions at high-redshift where remov- 
ing galaxies less luminous than Mb = —19.8 has a strong 
effect on the calculated tidal field. This is because some 
effort has been expended to get the redsh ifts and mag 



nitudes of galaxies along the line of sight (Morris et al 



1993), with a claimed completeness to wib — 19.0, so that 
one ought to be able to detect a galaxy of Mb = —17.0 at 
z = 0.065. 




Fig. 5. — The normalized mass function as extracted from the 
Schechter Luminosity function (using a = —1.2) and the TF relation 
(solid line), both in the B-band. The dotted line is the mass function 
corrected for sub-halos (see §5.4). 



5.5. Attributed Positional Errors 

For clouds associated with galaxies, there are two ma- 
jor sources of error in the assignment of a value for T at 
the cloud's location, galaxy redshift errors, and peculiar 
velocities, particularly of the cloud. We discuss these in 
order below. 

In this paper the main concern is with void clouds, which 
are only slightly effected by galaxy redshift errors since 
typical redshift errors, on order 20-30 kms , imply mod- 
est positional errors relative to void clouds (~ 260 to 400 
kpc). 

Cloud peculiar velocity errors stem from th e "triple- 
value ambiguity" noted by Fully fc Shaya (1984| ) in which 
a given LOS velocity may occur at three different dis- 
tances when the object is infalling toward a larger mass. If 
clouds are responding to the peculiar gravity of a galaxy, 
they may acquire maximum peculiar velocities not much 
greater than v p ss v c — V GK. (see Eq. 28) when the clouds 
are close to the galaxyf^ A peculiar velocity of this order, 

4 If peculiar velocities are a measure of the galaxy potential, then 
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under a strict Hubble flow, results in an attributed dis- 
tance, R a = V GJC/H a . The average over-density within 
this distance is, 



p{R a ) = 2p c 



GICR t 



2p c 



(42) 



For the Galaxy (assume v c = 220 kms 1 ), the truncation 
radius R t is scaled up from R% = 500 h^g kpc to 684 kpc. 
The over-density within radius R a is p ~ 0.46p cr it, where 
in this case, R a ~ 3.0 h^ 5 Mpc. The median position angle 
(the angle between the LOS and the vector from cloud to 
galaxy) for clouds distributed isotropically about a galaxy 
is 60°, leading to the conclusion that the average LOS 
velocity will be ~l/2 of v p . Therefore, low cloud velocities 
relative to a galaxy suggest either that the cloud is distant, 
or that its velocity has dissipated against the gaseous halo 
surrounding the galaxy. 

It is also possible that clouds which are actually distant 
(say, 2 to 3 Mpc from a large galaxy), but have a small 
impact parameter, may appear at very nearly the same 
velocity as the galaxy; a slight peculiar velocity may bal- 
ance the Hubble flow if the cloud is near the t urna round 
radius. Thus, clouds which are within R„ ~ V G7C /H of 
a galaxy may not have their positions accurately predicted 
by the redshift due to the peculiar velocity. Such a huge 
positional error is not likely for the average cloud in GDS 
since it would require a low-probability projection of the 
cloud at large radius R « 3 Mpc (where the cloud density 
is low) to a small projected radius R p . However, again, 
this kind of error cannot substantially effect the position 
of a typical void cloud, nor its attributed T. 

Small groups of galaxies may, by their virial velocities, 
introduce error into their attributed positions. We have 
seen that voids may be considered to begin at distances of 
« 2 Mpc from a single ~ C* galaxy (§2.1), but that the 
distance at which a cloud is a t the same degree of iso lation 
from a group will be larger. Bahcall fc Cen (1995 ) show 



that the one-dimensional velocity dispersion in groups and 
clusters can be approximated by 



o-los ~ (90 ± 20) N 1/2 kms" 



(43) 



where N g is the number of galaxies in the group. For 
N g = 10, we have o~los — 284 kms -1 , or a spread of 
Sr ~ 1. 9/175* Mpc on either side of the center. But the 
scaling relations (Eq. 29) show that the radius at which 
a given Twould occur is ~ IO 1 / 3 times farther than for a 
single galaxy. If we consider 2 Mpc as a plausible min- 
imum distance from a major galaxy for a void to begin 
(at 2 Mpc from an C* galaxy the tide Tr — 0.24), then 
this translates to a distance ~ 4.3 Mpc for the group. If 
the group is at the same redshift as the cloud, most of 
the group galaxies would have attributed radial positions 
within an angle of <^24° of the center of mass, as seen from 
the cloud's position, producing a distance variation of less 
than ~9%. The tidal field at a distance of 4.3 Mpc 
from this group would be about 86% of that calculated as- 
suming its galaxies were located at the group center. This 
error is significantly less than the error in an individual 

with the truncated isothermal halo, the peculiar velocity will be of 
order v p ~ 2 1 / 2 v c sfRt/R. 



galaxy. Further, since the real apparent magnitudes of the 
group should differ randomly, the error in the tidal field at 
the cloud should have a lower value from apparent mag- 
nitude error than a cloud subject to the influences of just 
a few. An error of this magnitude is insignificant in the 
present analysis. 

5.6. Tides and Cloud Statistics 

Sightlinc analysis results in a value for tidal field 7~(A) 
along each LOS (Fig. 4). The value of the specified tidal 
field limit separates the cloud population into two catalogs 
- those defined by T as an upper limit T u , and those with 
T as a lower limit T L . Histograms of clouds for various 
T are presented in Fig. 6, showing the former, C(T U ), 
as dotted lines, and the latter, C(T L ), as solid lines. In 
these plots, what is referred to as "void" is relative to the 
current choice of T u . Thus, "void" is relative to the stated 
tidal field limits. If a corresponding physical distinction 
between void and GDS clouds can be found, then "void" 
may take on a more absolute meaning. Note that the 
void histograms have a most probable Doppler parameter 
of b = 30 to 35 kms -1 over the whole range of T u . The 
high- velocity wing of the distribution suggests the presence 
in voids of more massive clouds, while the cutoff at low b 
gives the distribution a large skewness, with values ranging 
from 0.87 to 1.06 in voids 0.01 <; T Hm <; 0.1. For clouds 

defined by a lower bound T , as one increases T L , the flat 
distribution apparent for lower values of T L gives way to 
a peaked histogram and most probable Doppler parameter 
of b ~ 60 kms" 1 for T > 0.3. The skewness at T L = 0.3 
is -0.46. 




80 20 
b— value (km/s) 

Fig. 6. — Histograms of void (dotted) versus non-void (solid) 
clouds as a function of cloud Doppler parameter for various Tn m . 
"Void" is defined by the value of the current tidal field upper limit 
T u . Notice that the shape of the void cloud distribution remains 
much the same (though its size necessarily gets larger) as T u in- 
creases until Tu m ~ 0.3. 



In Fig. 7a, the fractional redshift coverage as a function 
of tidal field upper limit has been equated to a fractional 
volume which has T < T u . Notice that 90% of the uni- 
verse has T <; 0.16 (logT <; 0.8). Fig. lb shows the 
fractional redshift coverage as a function of EW 
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It is of interest to attempt to determine the point T v 
which actually separates void from GDS space. Figure 
8 shows a differential histogram extracted from Figs. 6, 
which shows the number and b-value of clouds added to 
the void population in the stated interval. Panel 8a shows 
a distribution quite similar to the overall void population; 
fig. 86 shows a relatively flat distribution; Fig. 8c shows 
a distribution with characteristics of the GDS population. 
This implies that the break between void and GDS envi- 
ronments is at 0.1 < T < 0.3, which brackets T = 0.16, 
the tide at which 90% of the universe has a lower tide. Re- 
call that 0.9 is the filling factor for the unshocked popula- 
tion (§1). Thus we provisionally assign log T v = — 0.8±0.2 
as the contour dividing void and GDS space for ISOT ha- 
los, bracketing a void filling factor 0.87 < /„ < 0.84. 
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Fig. 7. — Panel (a): The fraction of rcdshift space with T < Tu m 
plotted against \o%% im . This fraction can be considered to be a fair 
approximation of the volume filling factor for voids defined by T < 
Tn m . Note that tides logT < —2.5 may seriously underestimate 
the relative tidal field, as C galaxies farther than 7.5 Mpc may 
produce scalar tides of this order (see §5.2). Panel (b): The fraction 
of redshift space with the sensitivity function (in mA) less than 
Wum as a function of log W (in mA). 
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Fig. 8. — The differential histograms from Fig. 8; histograms of 
clouds added to C(T U ) for a T u = 0.01 to 0.1, b T u = 0.1 to 0.3, 
and c T u = 0.3 to 1.0. The clouds added in panel (a) appear to have 
roughly the same distribution as is found for T u = 0.01 in panel (a) 
of Fig. 8. Panel b is distinctly flatter indicating that clouds which 
are non-void may be beginning to enter the population. Panel c 
begins to show the distribution characteristic of non-void clouds, 
with a most likely Doppler parameter near 60 kms^ 1 . This trend 
shows that the tidal field level at which the character of added void 
clouds experiences rapid change is at 0.1 < T < 0.3. 



6. THE EQUIVALENT WIDTH DISTRIBUTION FUNCTIONS 

In §3 it was shown that the likely effect of tides on a 
population of uniform clouds is to disrupt those with den- 
sity p <; 2Tp crit (Eq. 19). However, clouds bound by dark 
matter are not expected to be of uniform density (see §2). 
Because of this, it is reasonable to suppose that the effect 



of tides on the distribution of absorber EWs may have its 
greatest observable effect on the slopes of the EWDF. 

We have seen that the redshift coverage at low sensitiv- 
ity is restricted. The EWDF removes this observational 
bias. The EWDF of clouds of a catalog C is calculated as 
the redshift frequency of clouds with rest-equivalent width 
greater than W, directly for each equivalent width limit 
W by noting the number of clouds of any equivalent width 
greater than W within the ranges of wavelengths for which 
the sensitivity functions have values less than W. This is 
then divided by the summed redshift range for which the 
sensitivity S (in ml) is smaller than W (Fig. 7b). Thus, 
the cumulative EWDF is, 



/(W) = 



<W(> W) _ Ew "ci(Wcj > W > S) 



dz 



(44) 



where it is read, "n c ; such that W c i > W > S" , and "Az 
such that W > S" . Note that while this distribution is 
ostensively cumulative, it is not impossible for the distri- 
bution to decline at low EW since there may be conditions 
in which the number of clouds in the redshift interval with 
very low sensitivity S is inordinately small. So while this 
distribution is cumulative, its calculation at each A is to 
some extent independent of its neighboring bins. 
Results are presented below. 
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Fig. 9. — Equivalent width distribution functions, expressed as 
the logarithm of the cumulative line density as a function of the log 
of the EW (in mA). Panel a, an analysis of the derived mean EWDF 
(heavy solid line), and that derived by PSSII (heavy dashed line), 
and the EWDF for C(T U ) = 0.01 for z < 0.036 (dotted). Thin lines 
are number-weighted linear fits to current data. Note the inflection 
in our mean EWDF (heavy line) at log(EW) = f .5, which is not seen 
in Penton et al. analysis, is continuous with the void EWDF. Panel 
b, EWDFs with Poisson errors. The mean distribution function 
(heavy solid line), and the lighter lines are the EWDFs for catalogs 
formed using Tu m = 0.1 as lower-bound (upper line) and upper- 
bound (lower line). 



6.1. Standard Processing 

Standard processing uses the ISOT model, accepts all 
4-er clouds, and attributes galaxy mass based on apparent 
magnitude tub, a flat- A model, and the TF relation perti- 
nent to the B-band. Fig. 9a shows standard processing for 
all data irrespective of tidal field (which we hereafter refer 
to as "mean"), shown as the heavy solid line. Also shown, 
closely associated, is the mean EWDF taken from PSSI 
(dashed line). The two are very close except at W ^ 80 
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mA; where the PSSI shows a broad bump, the present 
EWDF shows an inflection to a higher slope. This differ- 
ence is due to the different sensitivity functions employed 
in the present analysis. Also shown (the dotted line) is 
the EWDF for T < 0.1, where the data used excluded the 
range z > 0.036. Notice that its slope appears to be con- 
tinuous with that of the inflected part of the mean EWDF 
(solid line). This is what one would expect if the low-EW 
part of the mean EWDF is dominated by clouds in voids. 
Also shown (thin solid lines) are the number-weighted least 
squares linear fits to the EWDFs. 

In Fig. 96 we show the Poisson errors associated with 
points for the mean EWDF (heavy solid) and the EWDFs 
formed by the complementary catalogs defined by % im = 
0.01; GDS is above, void is below. The error bars are 
not independent, but simply reflect the number of clouds 
at that EW which are contributing to the EWDFs. The 
cutoffs evident in the EWDFs at high EW are probably 
due to the tendency of large EW clouds to have associated 
metal-line absorption, and so be rejected from the line 
lists. The cutoff in the C(T L ) EWDF (the upper line) 
at low EW is the result of the paucity of low EW clouds 
in higher tidal field environments. The error bars show 
that the two complementary populations are significantly 
distinguished. 

The processing of the EWDFs was done for a series of 
Tiim equally spaced in logT. Number- weighted {i.e., the 
number of clouds contributing to the EWDF at a given 
EW) linear fits to the equation, 

\og{dN/dz) = log/(W) = C + Slog(W/63), (45) 

are calculated, where C is the intercept \og(dAf jdz) at 
W = 10 18 ~ 63 mA (the approximate median logW) in 
the total cloud catalog, and S is the slope. Errors are 
derived from the residuals and the diagonal elements of 
the variance-covariance matrix without number weighting. 
These results are presented in Table 2. For comparison, 
the fits for the mean EWDF are also presented. Because 
of the break in the distribution apparent at W ~ 32 mA, 
the fit is calculated for W < 32 mA, all W, and W > 32 
mA, respectively. Notice that the fitted parameters for the 
first resemble fits for void EWDFs (T <; 0.1). 

The data from the lower two subdivisions of Table 2 
are plotted in Fig. 10, where panel a shows the trend 
of slope S with logT, and panel b shows the intercepts 
C. Slope fitting errors are large compared to intercept 
errors. The trends in Fig. 10 indicate that in higher 
tidal field environments the redshift frequency of absorp- 
tion systems is larger, and the population more uniform 
in their EWs, while for low tidal field environments the 
slopes are steeper (more negative), representing an abun- 
dance of low EW clouds. The trend in the void slope 
with logT shows signs of a relative discontinuity some- 
where between 0.1 <; T <; 0.3, the range within which we 
hypothesize that the "true" void leaves off and the GDS 
begins (§5.6). The same trend appears less obviously in 
the plot of the intercepts (panel b) . If this trend is real it 
would suggest that the deepest voids may have marginally 
more large EW clouds and marginally larger line densi- 
ties than in locations closer to void edges. Note also that 
choosing a larger lower tide limit (upper sub-panels) is 
choosing smaller volumes around galaxies, and results in 



Table 2 

Fitted slopes and intercepts as a function of T 

LIMITS 



Mean EWDF 


s 


os 


C 




W < 32 mA 


-1.6670 


0.0808 


1.723 


0.0345 


W > 12 mA 


-0.9665 


0.0285 


2.082 


0.0129 


W > 32 mA 


-0.8900 


0.0387 


1.949 


0.0166 


T u 


S 




C 


o~c 


0.010 


-1.4688 


0.0583 


1.7812 


0.0227 


0.030 


-1.5277 


0.0522 


1.7204 


0.0203 


0.100 


-1.5303 


0.0557 


1.7383 


0.0217 


0.300 


-1.2823 


0.0398 


1.8233 


0.0155 


1.000 


-1.1417 


0.0535 


1.8912 


0.0235 


3.000 


-1.0267 


0.0619 


1.9450 


0.0272 




S 


os 


C 




0.010 


-0.6690 


0.0904 


2.1087 


0.0402 


0.030 


-0.5974 


0.0939 


2.2487 


0.0417 


0.100 


-0.4845 


0.1207 


2.5602 


0.0528 


0.300 


-0.5234 
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significantly larger line densities (panel b). The trend in 
the line density in GDS is consistent with C oc T 43 . The 
tendency for large errors in attributed T to be introduced 
in the GDS via positional errors makes this value highly 
approximate. 

One may conclude that there is a significant variation in 
the EWDF between the upper- or lower-limit tidal fields. 
The systematic variation of EWDF slope with Tim is con- 
sistent with expectations of the effects of tidal fields (§3). 
The modestly slight variation in distribution functions for 
a range below T u w 0.16 indicates that the void cloud 
EWDF is fairly uniform over the ~ 90% of the universe 
occupied by voids. The increasing concentration of clouds 
around galaxies, as shown by the steady steep rise in the 
intercepts of catalogs C(T L ) as T L increases, appears con- 
sistent with infall of clouds into the potential wells of 
galaxies and groups of galaxies. 

6.2. Error analysis by the variational approach 

We need to assure ourselves of the robustness of these 
results. In previous sections a number of possible biases, 
uncertainties and sources of error were noted which could 
not be directly evaluated except by noting the effect of 
different treatment of the data on the derived EWDFs. 
The results of these comparisons are presented below. 

6.2.1. Completeness 

Recall (§5.4) that the effect of incompleteness is to put a 
skew on the depth of the luminosity function probed across 
a range of redshifts, and hence to under-estimate T at the 
high-end of the redshift range, relative to that at low red- 
shift. This effect is shown graphically by the dotted lines in 
Fig. 4, where an absolute magnitude upper-limit of —19.8 
was imposed on the galaxy catalogs. In this section, the 
data is separated by a redshift constraint, and the EWDFs 
are re-derived for the low-redshift part, and compared with 
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Fig. 10. — Standard model fits: The number-weighted linear fits 
to the slopes S (panel a) and intercepts C (panel b), terms defined 
in Eq. 45, are presented for various T limits of cloud catalogs (see 
text). The solid lines represent a linear fit to the slopes (panel a) 
and intercepts (panel b) as a function of the logarithm of the tide T. 
In panel a, the upper horizontal dotted line has a value of the slope 
of the mean EWDF for W > 32 mA (the location of the inflection 

in the mean EWDF), and the lower, W > 15 mA. Above that 
(filled circles) are the EWDF slopes for catalogs C(T L ) (non-void), 
and below the dotted line (open circles) are the slopes for C(T U ) 
(void clouds). In panel b, the horizontal dotted lines represent the 
intercept (logarithm of dj\f /dz) for the mean EWDF for W > 32 
mA (lower) and for VV > 15 mA (upper). Filled circles represent the 
intercepts of C(T L ) EWDFs, and open circles represent the EWDF 
intercepts for catalogs C(T U ). The short dotted lines in the void 
regions of panels a and b are linear fits to the points with T u < 0.16 
indicating a possible inflection of the trend of slope and intercept 
with T u . 

our standard. We make our division at z — 0.036, at which 
point a apparent magnitude limit of tub = 16 means that 
we are seeing all but galaxies with C < 0.74£* (v c ~ 145 
kms -1 ). Fig. 11 shows the results of that analysis in 
comparison to that of the standard processing for catalogs 
based on Tu m = 0.01 and 0.1. Figure 12 shows the lin- 
ear fits for all C(T) for z < 0.036. The lines are the fits 
to data at all redshifts (Fig. 10), placed for the conve- 
nience of a visual comparison of the changes. Note that 
the slopes of void clouds have steepened; at Tu m = 0.1, 
the slope of EWDFs rises from ~ —1.5 to ~ —1.7, pos- 



sibly a better estimate of the true EWDF slope in voids. 
The previous impression of an inversion of the trend of 
slopes for T u <; 0.16 (panels a and b, short dotted lines 
from Fig. 10) is now strengthened in significance (dashed 
lines). Though the number statistics are now weaker, it 
appears that we may have gained in some respects by dis- 
carding highcr-rcdshift data, for the feature at T ~ 0.16 
in Fig. 12 which was somewhat weak in the standard pro- 
cessing using all the data now appears stronger with the 
lower-redsthft data. This could be expected because fewer 
clouds from higher tidal field environments may contam- 
inate the low tide catalogs when higher redshift data are 
culled. 




Fig. 11. — Completeness test: The equivalent width distribu- 
tion functions (7ji m = 0.01 and 0.1, with (solid) and without (dot- 
ted), the redshift interval z > 0.036. The median EWDF is plotted 
(heavy solid line) for reference. The T u (void) EWDFs are evi- 
dently somewhat steeper with the redshift cutoff; for T u = 0.1, the 
slope steepened from -1.56 to -1.70. The T L EWDFs (upper pairs 
of lines) show a drop with the rejection of higher redshift data. This 
may be consistent with evolution effects. 



6.2.2. Effect of early-type galaxies 

Recall (§5.1) that in attributing mass to galaxies the 
TF relation was used to transform the absolute magni- 
tudes into galaxy masses. Early-type galaxies do not fol- 
low the TF relation, however, they do follow a quite sim- 
ilar relation, basically C cx er~ 4 . This is steeper than the 
TF relation, but the main effect, for our purposes, is that 
the mass-luminosity ratio is larger for early-type galaxies. 
This can be checked by alternatively using, or not using, an 
altered attributed mass for those galaxies listed in the CfA 
catalog as having negative de Vaucouleurs T-Types (those 
are generally lenticular or elliptical galaxies). The alter- 
ation involves multiplying the normally attributed mass 
by a factor of 3.0 when the T-Type is less than or equal 
to zero. Since ellipticals are more often found in clusters, 
we might expect that this would attribute more masses to 
large groups, or clusters of galaxies, inflating the T L con- 
tours somewhat, which would reduce the EWDF for the 
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Fig. 12. — Completeness test: Same as Fig. 10, except that 
slopes and intercepts are fit to catalogs C constrained by z < 0.036, 
and the solid lines are the previous fits found for catalogs with no 
redshift constraint for comparison. What appeared as a marginal 
trend for increase in slope for GDS has now disappeared, though 
the intercept relation remains relatively accurate (top sub-panels 
of panels a and b, respectively). Note that what was previously 
a marginal tendency for slopes and intercepts of void clouds with 
Tlim < 0.16 (short, dotted lines) to have a trend contrary to the 
general trend suggested by the linear fit, is now more significant 
(short, dashed lines). This provides evidence for a relatively sharp 
change of character of cloud distributions near Tn m ~ 0.16 with the 
truncated isothermal model. 



GDS at a given T L . There should be relatively little effect 
on the void EW DFs since clouds are thou ght to avoid high 
density regions ( Morris et al., 1991 , 1993 ) where early- type 



galaxies are concentrated. 

As can be seen in Fig. 13, there is slight reduction in 
the line density of the DFs. The reason for this is probably 
that the added mass simply increases the volume of space 
at a given tidal field limit, which then has a slightly lower 
EWDF, as noted in Fig. 106. 

6.2.3. Effects of halo type 

The NFW halo is significantly less massive than the 
truncated isothermal halo (§5.1, and Fig. 3), and therefore 
when halos of galaxies are interpreted as NFW, the main 
effect is that the tidal fields will be lower by a factor of ~ 5. 
Thus we can expect that the intercepts and slopes will be 




Fig. 13. — The T-Type test: The equivalent width distribution 
functions with (dotted lines), and without (solid), the T-Type cor- 
rections for early-type galaxies, which makes the attributed galaxy 
mass 3 times greater. Tide limits are Tn m = 0.01, 0.1, 1.0. The ten- 
dency is for EWDFs based on tidal field lower limits to be more in- 
fluenced in high tidal field, than low tidal field environments, though 
only marginally so. 



shifted to lower logT by a factor of w 0.7. These results 
are presented in the form of fitted slopes and intercepts of 
number-weighted linear fits in Fig. 14. 

It is difficult to judge from Figs. 10 and 14 alone which 
halo is most likely to be more correct. However, we saw 
signs in Fig. 8. on the basis of most probable b- values, 
that the transition from void to GDS clouds occurs in the 
range 0.1 < T < 0.3. A log(T) offset of -0.7 would place 
this transition at T w 0.03 (logT ~ -1.5). While nei- 
ther of these tides could not reasonably be responsible for 
truncating the clouds in the transition zone, at z = 1, a 
lookback time of 57% of the age of the universe, the ex- 
pected tidal field is ~ 8 times greater, which according to 
Eq. 19, constrains stable clouds to have densities greater 
than ~ 2.5, and 0A8p cr it for isothermal and NFW halos, 
respectively. However, it does not seem reasonable to ex- 
pect that the GDS and void clouds could be distinguished 
at a level as low as 0.5 p cr n level; clouds of such low density 
could not reasonably be expected to be detected with the 
current technology. Such a cloud in the Hubble flow would 
produce a column of only 10 12 cm~ 2 (W < 4 mA) over 
a velocity range of ~ 220 kms -1 . Hence, it appears that 
gas at a density of ~ 2 — bQbPcrit is being stripped at the 
void-GDS boundary, where the tidal fields by themselves, 
even in the last few Gyr, could not reasonably be expected 
to be effective, unless galaxy halos are significantly more 
massive than even the isothermal model used here. 

7. DISCUSSION 

Numerous 3-D cosmological simulations have suggested 
that voids have a low density, and would present very few 
absorbers above present-day obsesrving limits. However, 
careful separation of clouds into catalogs according to the 
ambient tidal field has resulted in the discovery of a void 
cloud population that is by no means insignificant. The 
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Fig. 14. — Halo-type test: Same as Fig. 10, except that slopes 
and intercepts are fit to catalogs C with tides generated by halos 
interpreted as NFW. The solid lines are the previous fits found for 
catalogs with isothermal halos. To first order the points are dis- 
placed by ~ 0.7 to the left (in logT) relative to the points in Fig. 
10 for the isothermal halo due to the smaller NFW halo mass. 



line density of void clouds reaches a value d/V jdz w 500 
for W > 10 12 ~ 15.8 mA. An order of magnitude calcu- 
lation may convey the meaning of this. Equation 4 shows 
that at z = 0, and r p — 50 kpc, the comoving density 
n ~ 16 Mpc -3 , and the mean free path between clouds 
I ~ 8 Mpc, or 2.4 A. An average cloud mass of 10 9 M Q (a 
reasonable mass for a just-detectable truncated isothermal 
halo qf circular velocity 
in voids: — However 



« c ^ 11 kins -1 ) implies fl m = 0.1 
an intci-cloud medium of unknown 

density must exist which may add substantially to the den- 
sity in voids. 

Supporting the contention that these are actually void 
clouds is the clear distinction between the slopes of EWDFs 
based on the complementary catalogs C distinguished by 
limiting tidal fields Tu m , and the difference in the most 
probable 6-values of their clouds. Under the truncated 
isothermal halo model, the transition from void to GDS 
environments occurs at logT k, 0.8 ± 0.2, or T ss 0.16. 
Not only do the slopes and intercepts of the void EWDFs 
undergo a strong change (Figs. 10, 12), but the differential 
histogram of cloud 6-values also shows the transition. At 



this tidal field limit, ~ 90% of the universe is void, and 
10% is GDS. 

That the change in cloud character noted above can be 
associated with a certain tidal field limit under a given as- 
sumed halo type supports the contention that tidal fields 
are a relevant quantity, associating cloud stability and 
cloud characteristics with void and GDS environments. 
It further supports the contention that clouds are self- 
gravitating to some extent, for the variation of EWDF 
with tidal field is consistent with the predicted effects on 
such clouds. In the case of void clouds, the lack of strong 
tides allows a distribution of clouds which is unconstrained 
to the current detection limits, allowing small, delicate 
structures to exist. However, their true nature cannot be 
known without successful modeling of them in cosmologi- 
cal simulations. In the case of GDS clouds, the tidal trun- 
cation which galaxies and groups of galaxies impose on 
clouds appears to have resulted in a paucity of low-EW 
clouds, and a general flattening of the EWDFs. On the 
other hand, if GDS environment is well-characterized by 
the filamentary structures seen in 3-D hydro simulations, 
and have a relatively high density and temperature, it is 
possible that much of the observed truncation is by ram 
pressure and shear forces as clouds fall into the GDS from 
voids. This would solve the problem (noted in §6.2.3) 
of tidal field levels at the location of the apparent tran- 
sition between void and GDS being too low with either 
halo model to easily reconcile with the probable density of 
gas being stripped from clouds moving from void to the 
GDS. Both effects have a common source in the gravita- 
tional potential of the mass concentration, and a similar 
result in the stripping of low-density gas, so that even if 
ram-pressure stripping is a major factor in the difference 
between void and GDS clouds, the tidal field is still an 
effective and valid way to separate void and GDS clouds. 

There is a hint of structure in voids - deep voids ap- 
pear to have marginally larger line densities, and shal- 
lower EWDF slopes; a hint of more massive clouds there. 
This may suggest some degree of hierarchical clustering 
has taken place. 

In comparing ISOT and NFW models, the most obvi- 
ous difference is that in mass (NFW is a factor of 4 to 5 
smaller). The NFW halo has a weaker ability to tidally 
truncate clouds, and to maintain the filling factor of 0.1 
for GDS, one requires that Tu m ~ 0.03 has in the past 
tidally truncated the GDS clouds (see §6.2.4). But this 
may be too small to reasonably have tidally truncated 
those clouds. In addition, the scaling of the ab sorption ra- 
dius with luminosity in the B-, and iif-bands ( |Chen et al 



2001) agree nicely with the scaling relations of ISOT 



Eq. 

sec 



29)7but do not fit wcll with the NFW halo ( Eq _ 4 

§§5.1.1,5.1.2). 

Though I consider these general findings to be secure, 
many improvements could be made in the quality and 
quantity of data. For instance, the depth and quality of the 
data in the CfA redshift catalog leaves much to be desired. 
However, when the galaxy catalogs from the Sloan Digital 
Sky Survey become available, much of these problems will 
be solved; both deeper (m s * = 17.65 mag), more consis- 
tent photometry, and generally superior redshifts will be 
available, so that tidal fields may be calculated to greater 
accuracy. 

The quantity of absorption systems, and the sensitivity 
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of the observations also leaves something to be desired. 
However, there will soon be published results of Space 
Telescope Imaging Spectrograph (STIS) observations of 
~ 120 new absorbers (M. Shull, priv. comm.), more than 
doubling the current number of absorbers. Somewhat far- 
ther in the future, the Cosmic Origins Spectrograph (COS) 
will provide significantly deeper spectrography of low red- 
shift AGN, and should greatly increase the number of HI 
absorbers available for study. 

While this analysis has displayed the nature of the dis- 
tinction between void, and GDS clouds, we know little 
about the clouds themselves. The EWDF is just a sam- 
pling of the clouds weighted by cloud cross-section at a 
given EW. We do not know their mass, their cross-section, 
or their density. We have a tentative cloud model - clouds 
of primordial gas whose expansion is restrained by dark 
matter halos - but do not know in any detail what par- 
ent population has given rise to these absorbers, nor what 
contribution they make to f2 m . Though analytical models 
may help to clarify what the current state of clouds might 
be, the accurate implementation of the cloud hypothesis 
adopted here requires a hydrodynamical study beginning 
at the point of reionization in order to assess the current 
extent and dynamical characteristics of void clouds under 
a variety of halo models. These results could be used to 
substantiate, or invalidate, those models. If the former, it 
would help determine the number and mass-distribution of 
clouds which are consistent with the observed void EWDF, 
and thus estimate their contribution to f2 m . This task re- 
mains for a paper now in preparation. 

8. CONCLUSIONS 

The large line density of clouds in voids, where simu- 
lations have predicted very few clouds at EWs accessible 
to even STIS observations (R. Dave, priv. comm. 2001), 
appears to require massive dark halos to restrain the dis- 
persal of gas needed to produce detectable absorbers. We 
have seen that such halos are consistent with the general 
predictions of hierarchical structure formation if small ha- 
los did not generally form Population III stars. We further 
find that the increase in line density of clouds with increas- 
ing tidal fields in GDS are consistent with continuing infall 
of clouds from void regions into the GDS. 

We make the following specific conclusions. 

1. Tidal fields are a convenient and effective parameter 
with which to separate clouds into the complemen- 
tary environments; "void", and "galaxy dominated" 
spaces. However, it is possible that the tidal field 
stripping of clouds is significantly supplemented in 
GDS by ram-pressure stripping by the diffuse, ambi- 
ent dissipated gas populating filamentary structures 
which appear in 3-D cosmological simulations. 

2. The EWDF is significantly steeper in voids than the 
median EWDF, suggesting that the void cloud pop- 
ulation has most of its mass at smaller scales. 

3. Most probable void cloud Doppler parameters are 
characteristically half that of GDS clouds. 

4. The line density in GDS is of order more than 10 
times larger for W > 63 mA, but for W <^ 35 mA, 
void clouds dominate the mean EWDF. 



5. The void cloud EWDF is essentially constant over 
- 90% of the universe, while the GDS EWDFs 
have line densities that increase dramatically with 
increases in T L . 

We have analyzed the void-cloud population and found 
significant differences between that and the picture pro- 
vided by 3-D hydrodynamic simulations. It is speculated 
that the lack of resolving power in the simulations may 
be partially responsible for this difference. However, it 
appears that even so, the sheer number density of clouds 
which must be in voids is startling, and may presage a 
shift in thought about Lya clouds. 
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port for this research, and helpful suggestions. I am also 
grateful for the many comments and helpful suggestions 
from Christopher F. McKee. I thank Steve Penton and J. 
Michael Shull for their help with the sensitivity functions 
and cloud identification issues. Finally, I wish to thank 
the anonymous referee for many helpful suggestions and 
comments. 

APPENDIX 

TRUSTWORTHINESS OF THE CfA REDSHIFT 
CATALOG 

A lot will depend on the reliability of the redshifts and 
magnitudes in the assembled catalogs. Since the CfA 
galaxy catalog is a compilation of data from many sources, 
its tabulated magnitudes are somewhat inhomogeneous. 
To check the validity of these magnitudes, the galaxy po- 
sitions are used to query the Automated Pl ate Scanner 
Facility at t he University of Minnesota {e.g., Pennington 
et al. (1993ffl The APS is based on digitized Palomar 



Observatory Sky Survey (DPOSS) plates. Not all galaxies 
were "found" by the APS. Reasons were often that the 
relevant plate was not mounted. Sometimes the galaxy 
was simply not found, though STScI digitized sky survey 
images^ clearly showed their presence (often these were 
edge-on spirals). In a very few cases, nothing was found. 
Roughly half of the CfA galaxies were matched by an APS 
galaxy. The apparent magnitudes are converted to abso- 
lute magnit udes, corrected for galactic extinction ( 5chlcgc| 
et al., 19*9^ ), and inclination (Fully et al., 1998). A scat- 
ter plot makes it clear that some of the magnitudes were 
anomalous. However, the overwhelming majority of galax- 
ies fit the relation, Maps = Mc/a + (0.65±0.5) mag, with 
a slope of 1.004 

The stated error of the AP S magnitudes for integr ated 
isophotal galaxy magnitudes ( Pennington et al., 1993| ) was 
estimated to be gaps — 0.3 to 0.5 mag, but a more recent 
estimate is 0.3 mag (R. Humphreys, priv. comm., 2001), 
reflecting the added sophistication of their processing. The 
observed dispersion is 0.5 mag in the APS-CfA plot, which 
is also the dispersions of APS and CfA magnitudes com- 
bined in quadrature. Thus, we estimate &cfA — 0.4 mag. 

The 0.65 mag offset indicates that the APS magnitudes 
are fainter. We shall see that much of this offset can be 

The APS databases are supported by the National Aeronautics 

University of Minnesota.. The 



ittp:/ /aps. umn.edu/ 



and Space Administration and the 
APS databases can be accessed at 

6 Thc Digitized Sky Surveys were produced at the 
Space Telescope Science Institute under the U. S. Gov- 
nrnmont grant NAH W-91 flfi and may be accessed at 



http : I J 'archive.stsci.edu / 1 cgi — bin/dss form 
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Fig. 1. — Panels (a) through (c) show luminosity function param- 
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eters. Panel (a) s^"-"r° Ml <"> a fame 
using Jl = QJZ5 QBjgjvton et al 2001] ; pole et al. 2001| ; |Mobashcr 
et al. 1993; |Tulry 198g| ). Panel (b) shows the faint-end slope a as a 
function of \^JJ~ Panel (c) shows the log of the normalization $* 
with h = 0.75. Panel (d) shows the slope (3 (see Eq. 9fi §5 1) pf tV. 



Pierce 



TullYiPisher relation as a function of spectral band ( Tully . 
2000j)T] 



understood in terms of the different spectral band in which 
the magnitudes are measured. The CfA magnitudes are 
generally B-band magnitudes with A e // ~ 4200 — 4400 
A. The APS magnitudes are POSS O-band with A e // ~ 
4000 A. A compilation of published studies of luminosity 



func tions (LF) in a range of spectral bands ( Blanton et al 
200l|; |Cole et al. 200l|; jVlobashor, Sharpies, fe Ellis 1993 



Tully 1983 ; |TuUy & Pierce 200q ) show how the defining 
parameters of the fit vary with the spectral band, revealing 
a tight relation between </>*, and M* (Fig. Ala, c) as a 
function of spectral band. Let us analyze the CfA versus 
the APS magnitudes in the light of this figure. Note first 
that as one chooses bluer bands, the galaxy luminosity is 
smaller at the "knee" of the luminosity function, and the 
normalization <f>* is larger. 

It is the Tully-Fisher relation ( [Tully et al., 1998| ) which 
allows us to use the LF data to get our answer. TF can 
be expressed, 

M band = M* band - 2.5 [{3 band {\og(2v c ) - log(2w )}] . 

Applying this once each to O and B magnitudes and sub- 
tracting them, we have the transform from O to B for a 
galaxy of circular velocity v c , 



B = 0-{0* -B*) + 2.5 {/? - log 



(1) 



where flband is the Tully-Fisher slope (see Fig. Aid), and 
Vq is a fiducial velocity of 158 kms" 1 . A linear fit of the 
Bt, R, and / bands (Fig. Aid) implies (3q ~ 2.83, and 
Pb ~ 2.87 to 2.89. M* and cf>* change relatively fast in 
these bands (Mq = -19.0 mag, c/)* ~ 0.029h 3 Mpc~ 3 and 



M* B = 19.3 - 19.7 mag, 4>* B ~ 0.026 - 0.023h 3 Mpc~ 3 ) for 
the 4200 A and 4400 A versions of B, respectively. 

Applying these values of M* to Eq. 1 shows a shift of 
+0.4 to 0.7 magnitudes for a galaxy with v c = 158 kins" 1 
going from B to O band magnitudes; O is "fainter". For 
galaxies with v c smaller than 158 kms -1 , there is only a 
slight skew in the sense that for small galaxies, the shift 
in attributed absolute magnitude is slightly larger. 

One can conclude that the CfA magnitudes are gen- 
erally consistent, and trustworthy within the constraints 
established above. 
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